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Abstract 

We study weak-field solutions having spherical symmetry in f{T) gravity; to this end, we solve 
the field equations for a non diagonal tetrad, starting from Lagrangian in the form f{T) = T +aT”, 
where a is a small constant, parameterizing the departure of the theory from GR. We show that 
the classical spherically symmetric solutions of GR, i.e. the Schwarzschild and Schwarzschild-de 
Sitter solutions, are perturbed by terms in the form oc and discuss the impact of these 

perturbations in observational tests. 
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I. INTRODUCTION 




Since the very discovery of the accelerated cosmic expansion [l|, ^ , and its conhrmation 
due to multiple observations 


[s-S], it has been customary to investigate theories that extend 
general relativity (GR), in order to get an agreement with the observations, without requiring 
the existence of dark entities. Hence, motivated by this instance, which recognizes that GR 
fails in describi^ gravity at large scales (consider also the old issue of the rotation curves of 
spiral galaxies |^), several theories have been proposed to generalize Einstein’s theory. Some 
of these modihed models of gravity are geometric extensions of GR, in other words they are 
based on a richer geometric structure, which is supposed to give the required ingredients to 
support the observations. 

As a prototype of this approach, one can consider the f{R) theories, where the gravita¬ 
tional Lagrangian depends on a function / of the curvature scalar R (see [?, 8| and references 
therein): when f{R) = R the action reduces to the usual Einstein-Hilbert action, and Ein¬ 
stein’s theory is obtained. Another example is given by the so-called f{T) theories, which 
have similarities and differences with respect to f{R). To begin with, they are based on 
teleparallel gravity (TEGR) js], where the gravitational interaction is determined by tor¬ 
sion, and the torsion scalar T appears in the Lagrangian instead of the curvature scalar. 
Furthermore, the underlying Riemann-Gartan space-time is endowed with the Weitzenbock 
connection (instead of the Levi-Givita connection), which is not commutative under the 
exchange of the lower indices, and has zero curvature while non-zero torsion. Actually, Ein¬ 
stein himself proposed such an alternative point of view on gravitation in terms of torsion 


and tetrads 


10| . In fact, in the TEGR picture the tetrads held is promoted to be the 


dynamical held instead of the metric tensor. In spite of these diherences, TEGR and GR 
have equivalent dynamics: in other words, every solution of GR is also solution of TEGR. 
However, when TEGR is generalized to f{T) by considering a gravitational Lagrangian that 
is a function of the torsion scalar, the equivalence breaks down [u, 12|. As a consequence 
f(T) theories can be considered potential candidates for explaining (on a purely geometric 
ground) the accelerated expansion of the universe, without requiring the existence of exotic 
cosmic huids (see e.g. jl^). 

While f{R) theories gives fourth order equations (at least in the metric formalism, while 


they are still second order in the Palatini approach, see e.g. 


8|), the /(T) held equations are 
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second order in the field derivatives since the torsion scalar is a function of the square of the 
first derivatives of the tetrads field. Furthermore, as for f{R) theories, the generalized TEGR 
disp lays additional degrees of freedom (whose physical nature is still under investigation 
14| ) related to the fact that the equations of motion are not invariant under local Lorentz 


transformations 


15| . In particular, this implies the existence of a preferential global reference 


frame defined by the autoparallel curyes of the manifold that solve the equations of motion. 
Consequently, even though the symmetry can help in choosing suitable coordinates to write 
the metric in a simple way, this does not give any hint on the form of the tetrad. As 


discussed in 


16j . a diagonal tetrad -that could in principle be a good working-ansatz for 


dealing with diagonal metrics- is not a good choice to properly parallelize the spacetime 
both in the context of non-fiat homogenous and isotropic cosmologies (Friedman-Lemaitre- 
Robertson-Walker universes) and in spherically symmetric space-times (Schwarzschild or 
Schwarzschild-de Sitter solutions). 

The case of the Schwarzschild solution and, more in general, of the spherically symmetric 
solutions in /(T) gravity, is particularly important, because these solutions, which describe 
the gravitational field of point-like sources, allow to test f{T) theories at scales different 
from the cosmological ones, e.g. in the solar system. Such class of solutions -both with 


diagonal and non-diagona’ 


years, see for instance 


17 


tetrads- haye been received much attention during the last few 
23| . Indeed, f(T) theories can be used to explain the cosmic 


acceleration and observations on large scales (e.g. via galaxy clustering and cosmic shear 
measurements 2^), but we must remember that since GR is in excellent agreement with 
solar system and binary pulsar observations 2^, every theory that aims at explaining the 


large scale dynamics of the Universe should reproduce GR in a suita 


e weak-field limit: the 


same holds true for /(T) theories. Recently, solar system data 2g l27| have been used to 
constrain /(T) theories; these results are based on the spherical symmetry solution found 
by lorio&Saridakis [261], who used a diagonal tetrad. In this paper we follow the approach 
described in to define a “good tetrad” in /(T) gravity - that is consistent with the 
equations of motion without constraining the functional form of the Lagrangian - and solve 
the field equations to obtain weak-field solutions with a power-law ansatz for an additive 
term to the TEGR Lagrangian, /(T) = T -|- aT”'. 

This paper is organized as follows: in Section [TTl we review the theoretical framework 
of f{T) gravity and write the field equations, whose solutions for spherically symmetric 
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space-times, in weak-field approximation, are given in Section IIIII Eventually, discussion 
and conclusions are in Sections m and El 

II. /(T) GRAVITY FIELD EQUATIONS 

We start by briefly discussing the f(T) gravity framework that leads to the held equations. 
To begin with, we point out that, in this scenario, the metric tensor can be viewed as a 
subsidiary held, and the vierbein held is the dynamical object whose components in a given 
coordinate basis e“ are related to the metric tensor by 

9,,u{x) = 7]abel{x)el{x) , ( 1 ) 

where rjab = diag(l, —1, —1, —1). Notice that latin indexes refer to the tangent space while 
greek indexes label coordinates on the manifold. Hence, the dynamics is obtained by the 
action^ 


S=-^Jf(T)ed‘x + SM, ( 2 ) 

where e = det e“ = ^y—det{g^^) and Sm is the action for the matter fields^. Here /(T) is a 
differentiable function of the torsion scalar T, which is defined as 

T = , (3) 

where the contorsion tensor is defined by 

SV = \ - V + V) + , (4) 

and the torsion tensor is 

T\.=4{9.el-d,et) . (5) 

Varying the action with respect to the vierbein e“(a;), one gets the held equations 

e-%{e + eJ>S/-d,{T)fTT + = 47rGe/-y, (6) 

^ We use units such as c = 1. 

^ Notice that many authors write the gravitational Lagrangian in the form T + f(T), thus denoting the 
deviation from GR by means of the function f{T): on the contrary, here f{T) is the whole Lagrangian. 
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where is the matter energy-momentum tensor and subscripts T denote differentiation 
with respect to T. 

We look for spherically symmetric solutions of the field equations, so we start from the 
metric 

ds^ = - r'^dO? , (7) 

where dVL^ = <76*^ -l-sin^ Odcf)^. Due to the lack of Local Lorentz Invariance, tetrads connected 
by local Lorentz transformations lead to the same metric - i.e. the same causal structure 
- but different equations of motions, thus physically inequivalent solutions. This means 
that, even in the spherically symmetric case, for which symmetry helps us in choosing the 
coordinates and the metric tensor in a simple form, it is quite complicated to do an ansatz 
for the tetrad held. In particular, for the symmetry and coordinates we are dealing with, it 
turns out to be a mistake to choose a diagonal form for e“^: it does not properly parallelize 
the static spherically symmetric geometry in the context of f{T) gravity. 

Then, with this caveat in mind, it is possible to derive the held equations for the non 
diagonal tetrad: 




/ e ^/2 0 0 0 \ 

0 sin 6 cos (j) sin 9 sin cj) cos 9 

0 —r cos 9 cos 0 —r cos 9 sin 0 r sin 9 

\ 0 rsin6'sin0 —rsin6'cos0 0 J 

following the approach described in 1^: in doing so, the functional form of the Lagrangian 
and the specihc form of the torsion scalar are not constrained a priori. We remark once more 
that diherent choices of tetras, while giving back the same metric represent diherent physical 
theories. In this work we are interested in a specihc tetras that does not lead to a constant 


torsion scalar. In such a theory Birkhoh theorem does hold, as shown in m while the most 
general vacuum solution is not of Scwharzschild De Sitter kind, as it happens with tetrads 
for which T = const. Then, it is worthwhile to investigate the features of the spherically 
symmetric solutions in this case for a generic Lagrangian, that nevertheless should admit 
the Schwarzschild solution when it reduces to the teleparallel equivalent of GR. The held 
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equations are 

- h'^-^ (2 - + r^e®<’'>T - 2rB'(r)) + 

+ ( 8 ) 

ffT) p-B{r) 

-^-h fr ^ 2 ~ + r^e^*'^^T — 2rp4'(r)) = Airp (9) 

Jt [—4 + — 2rA'(r) — 2rB'{r) + r^A'(r)^ — r‘^A'{r)B'{r) + 2r^A"(r)] + 

+2rfTTT' (2 + 2e-®('')/2 + rA'(r)) = 0 (10) 


where p, p are the energy density and pressure of the matter energy-momentum tensor, 
and the prime denotes differentiation with respect to the radial coordinate r. Moreover, the 
torsion scalar is 


T 




[1 + + rA\r)\ . 


( 11 ) 


III. WEAK-FIELD SOLUTIONS 


Exact solutions in vacuum (p = p = 0) and in presence of a cosmological constant 
(p = —p) of the above held equations are thoroughly discussed in ; here we are interested 
in weak-held solutions with non constant torsion scalar, i.e. T' = dT/dr ^ 0. 

Indeed, for actual physical situations such as in the solar system, the gravitational held 
is expected to be just a small perturbation of a hat background Minkowski spacetime. As a 
consequence, we write 

= 1 +A(r), = 1 + 5(r) (12) 

and conhne ourselves to linear perturbations. Moreover, we consider Lagrangians of suffi¬ 
cient generality, that we write in the form /(T) = T + where a is a small constant, 
parameterizing the departure of these theories from GR, and \n\ ^ 1. 


To begin with, we consider the case n = 2, that has been already analyzed in 26|]. From 
flSl lTOll we obtain the following solutions 

_C^ 

r 
Cl 


A{r) = -32 


a 


a 


B{r) =96 — + 


(13) 

(14) 
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where Ci is an integration constant. Then, on setting Ci = 2M, we get the weak-held limit 
of the Schwarzschild solution plus a correction due to a: 


ds^ = ( 1 


2M 


32 ^ ] df 


1 + — + 96 ^ ) dr^ 


r^dn^ 


Eventually, the torsion scalar turns out to be 


T(r) = - - 128- 


a 


.4' 


(16) 


(16) 


These results can be compared to those obtained in 


26[ |. where a Lagrangian in the form 


/(T) = T + aT^ was considered. While to lowest order approximation in both cases the 
perturbations are proportional to 1/r^, the numerical coefficients are different: this is not 
surprising, since the authors in 26|] solve different field equations. In particular they use 
a diagonal tetrad, which constrains the torsion scalar to be constant (see e.g. lfl| and 
references therein): however, the solution given in 2^ does not seem to have a constant 
torsion scalar, which makes it inconsistent. 

Likewise, if we look for solutions of the equations (j8 |) - ffT0l) with p = k, p = —/c, which 
corresponds to a cosmological constant, we obtain 


ds^ = 1 


2M W . 

-32— —Ar 

.2 3 


r 




dr — ( 1 + 


2M a W . 

-+ 96— + -Ar 

r 6 


dr^ — r^dVr 


(17) 


where we set /c = ^ and A is the cosmological constant. The torsion scalar is the same 
as flT^ . So, the weak-held limit of Schwarzschild - de Sitter solution is perturbed by terms 
that are proportional to a. 

The previous results can be generalized to the case of a Lagrangian in the form /(T) = 
T -I- ckT"-, and we get 


A{r) = - 


Cl 


„2-2n 


a- 


-)3n—1 


2n — 3 


-AC 


(18) 


„2-2n 1 

Bir) = — -I- a -(—3?7, -|- 1 -|- 2m?) + -AC (19) 

r 2n — 3 3 

In particular, if A = 0, we obtain vacuum solutions. Notice that, on setting Ci = 2M, we 
obtain a weak-held Schwarzschild - de Sitter solution perturbed by terms that are propor¬ 
tional to a and decay with a power of the radial coordinate, the specihc value depending on 
the power-law choosed in the Lagrangian. The torsion scalar is 

T(r) = ^ + 2ar-^^2^^ (n + 1) (20) 


7 












while the perturbation terms due to the deviation from GR are in the form 


23n—1 

where an = - - bn = 


Aa{r) = aanT 

23n—1 


2 -^", B^{r) = abnv^-^^ ( 21 ) 

(2n^ — 3n + l). A close inspection of the perturbation 


2n-3’ 2n-3 

terms reveals that they go to zero both when r —)■ cxo with n > 1 and when r —)■ 0 with 
n < 1. In the latter case, in order the keep the perturbative approach self-consistent, a 
maximum value of r must be defined to consider these terms as perturbations of the flat 
space-time background. 

We remark here that our linearized approach can be applied to arbitrary polynomial 
corrections to the torsion scalar: as a consequence, by writing an arbitrary function as a 
suitable power series, it is possible to evaluate its impact as a perturbation of the weak- 
field spherically symmetric solution in GR: the n-th term of the series gives a contribution 
proportional to 

It could be interesting to test the impact of the perturbations (|2l|l . To this end, we 
remember that it is possible to obtain the secular variations of the Keplerian orbital ele¬ 
ments due to general spherically symmetric perturbations of the GR solution, describing the 


gravitational field around a point-like mass, as one of us showed in 


For instance, the 


average over one orbital period of the secular precession of pericenter turns out to be: 
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Uba 


2n 


1 23^-1 (2 - 2n) (3 - 2n) Vl - ^ 

< cu >= -a -^^-^2-F 

4 {2n — 3) 




for n > 


for n < 


( 22 ) 


(23) 


In the above equations Ub, a, e are, respectively, the mean motion, the semi major axis and 
the eccentricity of the unperturbed orbit, while F is the hypergeo metric function. These 
relations can be used to constrain the parameters a,n, on the bases of the ephemerides 
data. 


IV. DISCUSSION 

It is useful to comment on the constraints one can infer for the parameters of our model 
from solar system data. But, before proceeding, it is important to emphasize a point about 
the tests of f{T) gravity. In theories with torsion, there is a sharp distinction between the test 








particles trajectories: autoparallels, or affine geodesics, are curves along which the velocity 
vector is transported parallel to itself, by the space-time connection; extremals, or metric 
geodesics, are curves of extremal space-time interval with respect to the space-time metric 


29l |. While in GR autoparallels and extremals curves do coincide and we can simply speak 


of geodesics, the same is not true when torsion is present. So, it is not 


irivial to dehne the 


actual trajectories of test particles. The results obtained by [2^ and |27|, together with the 
expressions fl2^ and fl2^ of the secular precession of the pericenter, strictly apply to the case 


of metric geodesics. According to us, this is a very important issue, that is often neglected 
in the literature pertaining to theories alternative to GR based to torsion: we will focus on 
this issue in a forthcoming publication 3^. In the same publication, we will constrain the 


parameters a and n, taking into account the recent da 
System provided by INPOPlOa {m, 32| and EPM2011 


:a of the ephemerides of the Solar 
Actually, perturbations in 


the form of power-law are present in different models of modihed gravity , and their impact 
on the Solar System dynamics has been analyzed, for instance, in [361-138 1. 

Bearing this in mind, it is possible to comment on our results and compare them to those 
already available in the literature pertaining to f{T) theories. In particular, because of the 
different choice of the tetrad, our solution, even in the case of a quadratic deformation of 
the TEGR Lagrangian, differs from the one found by lorio&Saridakis. Both corrections are 
proportional to l/r^, but they have different numerical coefficients. 

In particular, on substituting n = 2 in Eq. fl2^ . we obtain 

16a 


< u >= 


a'^rib (1 


(24) 


On the contrary, the corresponding expression obtained by lorio&Saridakis 2q is 

3a 

< >is= 


(25) 


a'^Ub (1 — e^) 

We see that they differ for a factor 16/3: the same happens to the constraints that can be 
obtained from our solution, by applying the approach described in |26[ and [27 1. 

In particular, lorio&Saridakis 2q derive constraints from the rate of change of perihelia 
of the hrst four inner planets, obtaining 


|A| < 6.1 X 10 


-42 -2 

m 


\a\ < 1.8 X 10^ m?. 
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Tighter results have been obtained in a subsequent paper [27|, where the authors consider 
upper bounds deriving from different phenomena: perihelion advance, light bending, gravi¬ 


tational time delay 39l-l42|. But the strongest constraints come from the perihelion advance. 


in particular from some supplementary advances constructed by considering that the effects 
due to the Sun’s quadrupole mass moment might represent possible unexplained parts of 
perihelion advance in GR Ml- This gives 


|A| < 1.8 X 10"^^ m-2 
|q;| < 1.2 X 10^ m^. 

The upper bound for the solution in Eqs. flTTD would be 3/16 smaller, that is |a| < 2.3 x 
10 

Eventually, we comment on the issue of the parameterized post-Newtonian formalism 
(PPN), in the framework of f{T) gravity. In order to test theories of gravity that give 
rise to detectable torsion effects in the Solar System, a theory-independent formalism that 
generalizes the PPN formalism when torsion is present was developed in 29| (see also 43|): 
starting from symmetry arguments, the metric and the connection around a massive body 
are perturbatively expressed in terms of dimensionless parameters related to the matter- 
energy content of the source, namely its mass and its angular momentum per unit mass. 
In doing so, the new parameters, which add to the original PPN ones, can be constrained 
by the experiments. Our results, however, cannot be directly described in this framework: 
an inspection of our solutions fl2Tll clearly shows that the perturbations are not related to 
the matter-energy content of the source, rather they depend on a, which parameterizes the 
departure of the f{T) theory from GR (see e.g. Eq. (4.2) in 43|). So, a new formalism is 
required to test the content of the Lagrangian by means of observations: in a sense, a can 
be considered a new post-Newtonian parameter of this formalism. 


V. CONCLUSIONS 

We studied spherically symmetric solutions in the weak-approximation of f{T) gravity. 
In particular, we started from a Lagrangian in the form /(T) = T with \n\ ^ 1, 

where a is a small constant which parameterizes the departure of these theories from GR, and 
solved the held equations using a non diagonal tetrad, showing that, to lowest approximation 
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order, the perturbations of the corresponding GR solutions (Schwarzchild or Schwarzchild- 
de Sitter) are in the form oc a These results can be used to evaluate the impact of 

the non linearity of the Lagrangian, for instance in the solar system. 

The case n = 2, corresponding to the Lagrangian /(T) = T + aT'^ has been already 
analyzed by [2^ and |27|, who used solar system observations to set constraints on the 
parameter a. It is important to point out that the latter results are based on the solution 
obtained by lorio&Saridakis 2^ , where a diagonal tetrad was used, which forces the torsion 
scalar to be constant: however, that solution does not seem to have a constant torsion scalar, 
which makes the consequent constraints not reliable. 

On the other hand, since because of the invariance properties of f{T) the choice of the 
tetrad field is crucial, we performed our calculations by using a more general non-diagonal 
tetrad, according to the prescriptions given in |^, and obtained a new solution of the /(T) 
quadratic model, for which the torsion scalar is not forced to be zero. 

We used the results already available in the literature to obtain the correct constraints 
from solar system data on the a parameter for the Lagrangian /(T) = T + even if we 
pointed out that the distinction between autoparallels, or affine geodesics, and extremals, 
or metric geodesics, is crucial in /(T) gravity, and deserves further investigation that we are 
going to carry out in forthcoming publications. 
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